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Abstract. We present first results on event-by-event elliptic flow fluctuations in 
nucleus-nucleus collisions corrected for effects of non-flow correlations where the 
magnitude of non-flow correlations has been independently measured in data. Over the 
measured range in centrality, we see large relative fluctuations of 25-50%. The results 
are consistent with predictions from both color glass condensate and Glauber type 
initial condition calculations of the event-by-event participant eccentricity fluctuations. 
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1. Introduction 



Elliptic flow is one of the key observables in understanding the early stages and the dynamics 
of heavy ion collisions. Large elliptic flow signals have been observed at the top RHIC energies 
in both Au+Au and Cu+Cu collisions [UllllS]. These results have been understood in terms of 
the initial anisotropy in the collision region, best described by the participant eccentricity [3ll3], 
being preserved via an early thermalisation after the collision and the hydrodynamic expansion 
of the system with very little viscosity [S]. Initial measurement of V2 fluctuations appear to 
confirm this picture where the event-by-event fluctuations in the initial shape anisotropy, 
quantified by the participant eccentricity are translated to event-by-event fluctuations in the 
V2 signal f6|. 

Particle correlations other than flow (non-flow correlations) such as HBT, resonance 
decays and jets can resemble correlations due to elliptic flow and can have various effects on 
different flow measurements. In particular, non-flow correlations can broaden the observed 
event-by-event V2 distribution and could be mistaken for V2 fluctuations. Quantitative 
understanding of the effect of non-flow correlations on the observed V2 fluctuations requires 
a measurement where the flow and non-flow correlations are completely disentangled. In this 
paper, we introduce a new method of measuring the non-flow correlations which crucially 
relies on the large pseudo-rapidity coverage of the PHOBOS octagon multiplicity detector. 
We present results on non-flow correlation strength and elliptic flow fluctuations corrected for 
non-flow correlations in Au-|-Au collisions at y/s~^ = 200 GeV. 



2. Measurement of non-flow correlations 



We aim to measure the second Fourier coefficient of two-particle angular correlations and 
separate contributions from flow and non-flow in this signal. Flow and non-flow contributions 
can be separated with a detailed study of the rj and Ar] dependence of the Acp correlation 
function. To this end we define the two particle correlation function between particles from 
two r] windows centered at r]i and r]2: 

i?„(A(/.,r?i,r/2 = — — r-1, (1) 

B{A(f),r]i,r]2) 

where F is the foreground distribution determined by taking hit pairs from the same event 
and B is the background distribution constructed by randomly selecting particles from two 
different events with similar vertex position and centrality. In practice, the correlation function 
is calculated using the procedure described in Ref. |^ , correcting for the effect of the particles 
produced due to the interactions with the detector. 

It can be shown that if the only correlations between particles are due to elliptic flow, the 
correlation function will take the form: i?„(A(/>, ?7i, 772) = 2^2(^1) x V2{ri2) cos(2A(/>). In general, 
both flow and non-flow effects are present and so we denote the second Fourier coefficient of 
Rn{A(p,r]i,r]2) as 2^2(771, ?/2)- The seperate flow and non-flow contributions are denoted f2(^) 
and 5{r]i,r]2) respectively, such that: 

V2{vi,m) =v2{m) X V2{r]2) + Km,V2) (2) 

In Fig. [l| viim ,772) is shown for a selected centrahty range in Au-|-Au coUisions at 
^Sj^^ = 200 GeV. The contribution of two different sources can be observed: a flow plateau, 
V2i'ni) X ^^2(^2)) which is seperable in rji and r/2 and a non-flow ridge along the diagonal rji = rj2. 

At large pseudo-rapidity separations, i.e. Ar/ = \r]i — t]2\ > 2, we expect the non-flow 
component 5{r]i, 772) to be small and estimate it by a comparison of the Fourier expansion of 
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Figure 1: Second 
Fourier coefficient of 
the correlation function 
i?„(A0, 7/1,772) as a 
function of rji and 772 
for the 40-45% central 
Au+Au collisions at 
^ = 200 GeV. 



Figure 2: Circles: To- 
tal measured fluctuations 
including flow fluctuations 
and effects of non-flow cor- 
relations [6j. Squares: Ex- 
pected measured fluctua- 
tions for the observed non- 
flow correlation signal. 



Figure 3: Circles: Results 
corrected for effects of non- 
flow correlations. Thick 
black line ands squares 
show cr{epart)/{epart) cal- 
culated in CGC [9j and 
Glauber MC [fij models, 
respectively. 



the A(j) correlations between data and HIJING. From these studies we estimate roughly that 
the non-flow in data in this region is in the range to 3.2 times the non-flow magnitude in 
HIJING. The uncertainty in the estimation of S{rji,rj2) for At] > 2 is the dominant source of 
error in the final magnitude of non-flow correlations. Further studies are planned to improve 
the understanding of the non-flow magnitude in this region. 

Performing a separable fit of V2{r]i) x V2{r]2) to fK??!, ^2) — ^(vi^'H^) in the region Ai] > 2, 
where r/2) is measured in data and 5{r]i,r]2) is taken as 1.6 it 1.6 times the magnitude in 

HIJING, we can measure V2{ti). The fit in the selected A-q region can be used to extract the 
correlation magnitude due to fiow, V2{r]i) x ^2(^/2), in the whole pseudorapidity accceptance 
\r]i\ < 3 and \r]2\ < 3. Subtracting the correlations due to fiow from f2(7?i,f/2) we can calculate 
(5 (771,772) in all regions of 771 and 772- 



3. Elliptic flow fluctuations corrected for non-flow correlations 

If there are no true fiow fiuctuations, the measurement of event- by-event V2 fiuctuations will 
yield an RMS width of c"(^^2)non-fiow ~ \fW)T^^ where (J) is the average cos(2A(/)) due to non- 
fiow correlations over all particle pairs [8]. Since the second Fourier coefficient of Rn{A(f), 771, 772) 
is equal to 2(cos(2A(/))), the t^l defined above equals (cos(2A(/))). Therefore we can calculate 
{5) by averaging 5(7/1,772), the non-fiow component of V2 for particles in the two windows 
centered at r/i and r/2. Averaging over all particle pairs gives: 

The expected relative fiuctuations due to non-fiow, (^{v2)-^q^_^q-^/ (^2) is shown in comparison 
to the the total measured relative V2 fiuctuations in Fig. [2} 

When we study the response of the V2 fiuctuations analysis on event samples with fiow 
fiuctuations and non-fiow correlations, we find the following empirical dependence: 



^(^2)total = ^(^'2)fiow + ^(^'2)non-fiow ^ ^xp (^-^ 



(^2) non-fiow 
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where c''('^2)total ^^^^^ to the measured V2 fluctuations and 0'(v2)fjQ^ refers to the true flow 
fluctuations. Subtracting the expected fluctuations due to non-flow correlations, {5) /2, from 
the measured total V2 fluctuations via Eq. |4j we extract the corrected elliptic How fluctuations, 
shown in Fig. |3j 

Also shown in Fig. [s] are o"ep^rt/(^part) at fixed values of Npart obtained in MC Glauber [6] 
and color glass condensate(CGC) [9] calculations. The 90% confidence level systematic errors 
for MC Glauber calculations (shown as a band in Fig. [s]) are estimated by varying Glauber 
parameters as discussed in Ref. [3]. The relative flow fluctuations are in agreement with 
participant eccentricity predictions calculated both with MC Glauber and CGC type initial 
conditions over the full centrality range under study. The observed agreement suggests that 
the fluctuations of elliptic flow primarily reflect fluctuations in the initial state geometry and 
are not affected strongly by the latter stages of the collision. 

These results are therefore qualitatively consistent with a picture of the collision 
process in which the shape of the initial stage geometry follows the predictions of the 
participant eccentricity model and where the initial geometry is translated into the final state 
azimuthal particle distribution in a hydrodynamic expansion, leading to an event-by-event 
proportionality between the observed elliptic How and the initial eccentricity. The results, 
however, do not allow us to distinguish between Glauber and color glass condensate type 
initial conditions. 



4. Conclusion 

In summary, we have presented the first measurement of elliptic flow fluctuations corrected for 
non-flow effects extracted independently from data in Au-|-Au collisions at ^/s~^ = 200 GeV. 
A significant non-flow correlation strength is observed. However, even after correcting for 
these correlations, the elliptic flow fluctuation signal is seen to be large with a relative 
magnitude of 25-50%. We show that the magnitude and centrality dependence of these 
fluctuations are in agreement with predictions for fluctuations of the initial shape of the 
collision region quantified by the participant eccentricity calculated with either Glauber or 
color glass condensate type initial conditions. These results support conclusions from previous 
studies on the importance of geometric fluctuations of the initial collision region postulated to 
relate elliptic How measurements in the Cu-|-Cu and Au-|-Au systems |3]. The initial geometry 
seems to drive the hydrodynamic evolution of the system, not only on average, but event-by- 
event. 
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